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Abstract: We numerically study the dynamics of false vacuum bubbles which are inside 
an almost flat background; we assumed spherical symmetry and the size of the bubble is 
smaller than the size of the background horizon. According to the thin shell approximation 
and the null energy condition, if the bubble is outside of a Schwarzschild black hole, unless 
we assume Farhi-Guth-Guven tunneling, expanding and inflating solutions are impossible. 
In this paper, we extend our method to beyond the thin shell approximation: we include the 
dynamics of fields and assume that the transition layer between a true vacuum and a false 
vacuum has non-zero thickness. If a shell has sufficiently low energy, as expected from the 
thin shell approximation, it collapses (Type 1). However, if the shell has sufficiently large 
energy, it tends to expand. Here, via the field dynamics, field values of inside of the shell 
slowly roll down to the true vacuum and hence the shell does not inflate (Type 2). If we add 
sufficient exotic matters to regularize the curvature near the shell, inflation may be possible 
without assuming Farhi-Guth-Guven tunneling. In this case, a wormhole is dynamically 
generated around the shell (Type 3). By tuning our simulation parameters, we could find 
transitions between Type 1 and Type 2, as well as between Type 2 and Type 3. Between 
Type 2 and Type 3, we could find another class of solutions (Type 4). Finally, we discuss 
the generation of a bubble universe and the violation of unitarity. We conclude that the 
existence of a certain combination of exotic matter fields violates unitarity. 



Keywords: Black Holes, Classical Theories of Gravity. 



* jakobidetsortehul@gmail.com 
^eastone83@gmail.com 
■'■innocent© muon . kaist . ac . kr 



Contents 



i 



Introduction 

Basic results of the thin shell approximation 
|2.1 | Thin shell approximation 
^ Buildability of initial states and Farhi-Guth-Guven tunneling 
Conclusions of the thin shell approximation 



|3|. Beyond the thin shell approximation 
|3.1| Setup 

^2] The initial value problem and integration schemes 



3.3 Free parameters 



Causal structures and physical issues 

ITI Type 1: a collapsing shell 

|4.1.1| Collapsing shell solutions 
|4.1.2| Transition from stable to unstable field values 
Type 2: an expanding shell with unstable field values 
Expanding shells with unstable field values 
Dependence on potentials 
Stability analysis 



4.2 



4.2.1 



4.3 



Type 3: an inflating shell 
4.3.1 Physics of exotic matters 
The iV-shell bubble 
Simulations of A^-shell bubbles 



4.3.3 



4.4 



Type 4: a bursting shell 

Transition from Type 3 to Type 2: a bursting shell solution 



4.4.1 










11 

12 
12 
13 
17 
17 
17 
19 

3 
3 

n 



|5|. Discussion 30 
|5.1| Beyond the thin shell approximation ^ 
^ Generation of a bubble universe: discussion on the information loss problem 33 



|A] . Convergence and consistency tests 



A.l Type 1 



A.2 Type 2 



35 

36 

m 



A. 3 Type 3 and Type 4 



- 1 - 



1. Introduction 



Our Universe experiences an accelerated expansion; also, our Universe seems to have had 
a period of exponential expansion, so called inflation Q. The simplest explanation of 
these phenomena are to assume that our Universe is or was in a vacuum with non-zero 
vacuum energy Q, as the vacuum energy effectively gives a cosmological constant. After 
the discovery of landscape |^], it becomes natural to assume that there are a lot of different 
vacua. 

Once a complex potential structure is allowed, it is inevitable to include some kind of 
quantum tunneling from one vacuum to other vacuum. Some authors have studied these 
phenomena Q . For simplicity, we assume that the background is a positive false vacuum 
with a small cosmological constant or an almost flat background. This de Sitter space 
violates the energy conservation; thus any kinds of tunneling can be allowed in principle, 
i.e. tunneling from low to high vacuum is possible and of course, vice versa ||5|. 

If a true vacuum bubble is generated in the de Sitter background, then the true vacuum 
bubble will expand and dominate all over the background 0]. Here, we need to observe 
the causal structure of de Sitter (inside)-de Sitter(outside) combination Q. However, what 
will happen if a false vacuum bubble is generated? If the false vacuum bubble collapses, 
we need to observe the de Sitter-Schwarzschild-de Sitter combination 0[P|l^, where the 
Schwarzschild implies a Schwarzschild black hole due to the collapse of the transition region. 
On the other hand, if the size of a false vacuum bubble becomes on the order of the horizon 
size of the background ^] again we need to observe the de Sitter-de Sitter combination. 

The most difficult situation is when the size of a false vacuum bubble becomes less than 
the horizon size of the background but greater than the horizon size of the false vacuum 
\W. In this case, the false vacuum bubble has to inflate in physical coordinates but the 
outside observer will only see the Schwarzschild structure. Then, if we assume the null 
energy condition, this structure seems to be a kind of de Sitter-Schwarzschild-de Sitter 
space, where the Schwarzschild means the Schwarzschild wormhole of course, 

it is a mathematically allowed solution. However, when the null energy condition is violated 
1 12], its exact causal structure is not well-known, but maybe a dynamical generation of a 
wormhole will be accompanied. 

Next natural question is the meaning of the inside de Sitter space. It implies the 
generation of a bubble universe which is separated from our Universe, while an outside 
observer sees a black hole 0. Is it possible to happen in a laboratory? If we assume the 
null energy condition and global hyperbolicity, in a general relativistic sense, it seems to be 
impossible since the initial condition needs a kind of singularity from a singularity theorem 
P^ . This kind of bubble is known as an unbuildable bubble In a semi-classical sense, 
it seems to be possible via tunneling from a buildable bubble to an unbuildable bubble, i.e. 



tunneling from the outside to the inside of a Schwarzschild wormhole |14|||8|[pX|]. 

However, if the background is the anti de Sitter space and if one assumes AdS/CFT 
|16], false vacuum bubbles are expected to evolve by a unitary way. Therefore, one may 
guess that tunneling from a buildable state to an unbuildable state should be excluded in 
the anti de Sitter background |15||10|. Also, it is reasonable to apply this principle to a 
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background de Sitter space. 

Up to now, these results were based on the thin shell approximation. According to 
the thin shell approximation, the transition region needs a kind of energy shell to satisfy 
the Einstein equations. However, in a real situation, the transition layer will have a non- 
zero thickness and the field values of inside of the shell will have non-trivial dynamics. Of 
course, it is very difficult to solve such dynamics by hand, and hence one needs a numerical 
approach |17][18||19|. If we extend our calculations beyond the thin shell approximation. 



then we can describe not only the geometry (metric and shell), but also the field dynamics 
which was ignored in the thin shell approximation. 

In this paper, we prepare a false vacuum bubble inside of an almost flat background. 
According to the thin shell approximation, as one tunes the initial parameters, basically two 
behaviors are expected: namely the collapse or expansion of a shell 0]. We could reproduce 
the collapsing solution easily. However, one interesting issue is whether the expansion of 
a false vacuum bubble is possible or not. According to the thin shell approximation with 
the null energy condition, this seemed to be impossible unless assuming tunneling from a 
buildable state to an unbuildable state [0] [^] . 

Here, in this paper, we remark two important points: 

• First, an expanding bubble solution which contains the inflating region is difficult to 
obtain not by the reason of geometry, but by reasons of field dynamics. Field values 
of a false vacuum bubble is unstable as the shell expands. 

• Second, general relativity does not exclude a generation of an inflating bubble if one 
assumes exotic matter fields. We show that a wormhole is dynamically generated to 
induce a bubble universe along the shell. It does not necessarily require tunneling 
from the outside to the inside of a Schwarzschild wormhole. 

Since a generation of an inflating bubble implies the violation of unitarity, the authors 
suspect that some holographic arguments on the unitary evolution seem to have potential 
dangerous; or some holographic arguments restricts our assumptions on the initial state of 
a bubble. 

This paper is structured as follows. In Section ^, we study previous results of the thin 
shell approximation. In Section |^, we introduce numerical setup to extend beyond the thin 
shell approximation and introduce simulation parameters. In Section ^, we observe and 
classify solutions of false vacuum bubbles and discuss interesting physical issues. Finally, in 
Section ^, we comment our contributions from studies beyond the thin shell approximation 
and discuss the unitarity issue with the generation of a bubble universe. 

2. Basic results of the thin shell approximation 
2.1 Thin shell approximation 

We assume spherical symmetry and observe the dynamics of a false vacuum bubble inside 
of the true vacuum background. Traditionally, some authors studied this problem by using 
the thin shell approximation |0] 0] ||8| ||T^ . The thin shell approximation is based on 
the following two assumptions: 
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1. The inside false vacuum region is governed by the de Sitter metric and the outside 
true vacuum region is governed by the Schwarzschild metric. 

2. Between the two regions, there is a thin mass shell which has a proper surface tension. 
The metric ansatz for the inside and the outside regions, respectively, are 

dsf = -Mr)dtf + -j^dr^ + r^dn^, (2.1) 
Mr) 



and 

dsl = -fo{r)dtl + -j^dr^ + rW, (2.2) 
fo[r) 

where /i(r) = 1 — r^//^ and fo{i") = 1 — 2m/r. 
The metric ansatz for transition region is 

'^Ssheii = -dr^ + R{r)^dn\ (2.3) 
where r = R{t) holds. We get the equation of motion for the shell 



i?2 + f,(R) - ^i?2 + ;^(^) = (2.4) 

where k = Attct and a is the surface tension of the transition region. This can be reduced 
to the following form: 

+ VediR) = 0, (2.5) 

where 



Ves{r) = fo{r) . (2.6) 



However, to maintain the information of the sign of each roots, we need to compare the 
extrinsic curvature for the outside and the inside of the shell. The extrinsic curvatures are 
defined as follows: 

ft.m^^g)±^ = ±y]^I7^. ,2.7) 

and 

Po = f^(^^-fo(RJ--'^' = ±4¥7m). (2.8) 

Now, to satisfy the Einstein equations, 

(3i- l3o = kR (2.9) 

should hold. The sign of the extrinsic curvatures in i? — > or i? ^ oo limit intuitively tell 
us the asymptotic direction of the (expanding or collapsing) shell [ p!o|] . 
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singularity 




de Sitter Schwarzschild 



Figure 1: Solutions of the thin shcU approximation (symmetric cases). Since f3i is always positive 
in i? limit, dSn is disallowed; (3o is always positive in i? — ^ limit, Schc is disallowed; Po is 
always negative in i? — > cx3 limit, ScIie is disallowed. 



singularity 




de Sitter Schwarzschild 

Figure 2: Solutions of the thin shell approximation (asymmetric cases). 



In general, the effective potential T4fj is a convex function for a time-like shell [7| |lO| ; 
therefore, it allows a collapsing solution or an expanding solution. Then there are basically 
5 possibilities: (a) from expanding to collapsing, (b) from collapsing to expanding, (c) 
from collapsing to collapsing, (d) from expanding to expanding, and (e) a static solution 
in an unstable equilibrium, (a) and (b) are symmetric solutions, whereas (c) and (d) are 
asymmetric solutions. For simplicity, we omit the unstable equilibrium case (e). 

Firstly, let us classify symmetric solutions. The left diagram of Figure |l| is for the de 
Sitter space, and the right diagram is for the Schwarzschild space. For a collapsing case, 
dSA or dSD are possible; and Sche, Schc, or Scho are possible. Also, for an expanding 
case, dSe or dSc are possible; and SchA or Schs are possible. However, according to 
the behavior of the extrinsic curvatures in ii —> or i? —> oo limit, we can remove the 
solutions of dSn, Schc, and Schs- Therefore, there are 4 possible solutions: dSA — Sche, 
dSA — Scho, dSe — SchA, dSc — SchA. The case dSA — Sche is a collapsing bubble solution, 
where the collapsing shell is inside of a Schwarzschild wormhole. The case dSA — Scho 
is a collapsing bubble solution, where the collapsing shell induces a Schwarzschild black 
hole. The case dSs — SchA is an expanding bubble solution, where the shell expands inside 
of a Schwarzschild wormhole, and the shell becomes greater than the horizon size of the 
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singularity 




de Sitter Schwarzschild 



Figure 3: Farhi-Guth-Guven tunneling. 



inside de Sitter space. The case dSc — SchA is an expanding bubble solution, where the 
shell expands inside of a Schwarzschild wormhole, and the shell expands outside of the 
cosmological horizon for the r = observer. 

Secondly, let us classify asymmetric solutions (Figure |2|). The most interesting case 
is the creation of a bubble universe. In this case, we need to consider from expanding 
to expanding solution. Here, dSs, dSp, and Schp are allowed; thus giving us the case of 
dSE — Schp and dSp — Schp as allowed transition solutions. We can interpret these as 
expanding solutions which begin from a singularity |M]|13]. 



2.2 Buildability of initial states and Farhi-Guth-Guven tunneling 

If we want to build a bubble universe in a laboratory, the initial state should begin from 
the inside of the horizon of the de Sitter space, and the outside of the horizon of the 
Schwarzschild black hole or the Schwarzschild wormhole. And the final state should end 
at the outside of the horizon of the de Sitter space and the inside of the Schwarzschild 
wormhole. Then the only reasonable solution is the unbound solution dSs — Schp or 
dSp — Schp. 

However, this is not the final answer to the generation of a bubble universe. According 



to Farhi and Guth [13|, whenever the null energy condition and global hyperbolicity hold, 
if a shell becomes greater than the horizon size of the inside de Sitter space, the horizon 
becomes a kind of anti-trapped surface, and the bubble should begin from an initial sin- 
gularity. Therefore, if a solution ends with a bubble universe, its initial state should be a 
singular state. We call this an unbuildable state and the opposite case is called a buildable 
state ilOi. 



To overcome this problem and generate a bubble universe by using a constructible way. 



we need to include tunneling |14|. The tunneling is to paste a buildable solution dSA — Scho 
and an unbuildable solution dSe — SchA |§. The former does not hold the conditions of 
the singularity theorem, and thus, one may assume that its initial state is buildable. Some 
authors have calculated the probability using some approximations of quantum gravity, and 
found consistent results This is known as Farhi-Guth-Guven tunneling (Figure ^). 



-6- 



2.3 Conclusions of the thin shell approximation 

The followings are conclusions of the thin shell approximation with the null energy condi- 
tion: 

1. A shell will either collapse or expand; if it has sufficient energy, it can expand to an 
unbounded size. 

2. A shell cannot expand except if it is in a Schwarzschild wormhole. 

3. When a false vacuum bubble is generated in an almost flat background, if one wants 
to make a bubble universe, it should tunnel into the Schwarzschild wormhole. 

One interesting question is this: is it possible to prepare any field configurations, thin 
or thick shell, so that the bubble is able to expand forever? Let us assume an initial 
bubble which is inside of the de Sitter horizon and outside of the Schwarzschild black hole. 
Then, the only allowed solution from the thin shell approximation is dSA — Schn, i.e. a 
collapsing bubble solution P], and hence it cannot expand forever according to the thin 



shell approximation [^|1C]. Now the next natural step is to move to beyond the thin shell 
approximation. 



3. Beyond the thin shell approximation 

In this paper, by going beyond the thin shell approximation, we observe some field config- 
urations where the bubble which is outside of the Schwarzschild radius can expand forever. 
Then, what is the difference between the thin shell approximation and beyond the thin 
shell approximation? 

We specify the following key assumptions for beyond the thin shell approximation: 

1. The field of the inside region is in a false vacuum and the field of outside is in a true 
vacuum; we do not assume special metric structures and we consider whole dynamics 
of metric and fields. 

2. The transition region has a non-zero thickness. 

One possible interpretation is that, from the first assumption, the field value of the 
inside region is not static and not stable; then, the inside is no longer an exact de Sitter 
space. The other possibility is that, from the second assumption, the thick transition 
region induces a wormhole in a dynamical way. In this paper, we will demonstrate that 
(1) if we do not violate the null energy condition, the first possibility happens and (2) if 
we do violate the null energy condition, both possibilities happen. 

Now we discuss our model and setup. 

3.1 Setup 

We describe a Lagrangian with a scalar field $ with a potential V{^) pl|] p^: 

>C = -^';a^;65"' - 2F($). (3.1) 
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From this Lagrangian we can derive the equations of motion for the scalar field: 

^;af>5"'-^'W = 0. (3.2) 

Also, the energy-momentum tensor becomes 

Tab = ^;a«';6 " \g ab{^ .c^ ■,d9"' + 2V{^)). (3.3) 

Now, we will describe our numerical setup. We start from the double-null coordinates 
(our convention is [u,v^9,ip\), 

ds"^ = -a^{u, v)dudv + r'^{u, v)dn^, (3.4) 

assuming spherical symmetry. Here, u is the in-going null direction and v is the out-going 
null direction. 

We define main functions as follows (we follow the numerical approach of previous 
authors |23]|17]|18]|jl^.): the metric function a, the area function r, and the massless 
scalar field S = ^/Att^. Also, we use some conventions: d = a^^/ct, h = a^u/ct, f = f,u, 
g = r „, W = 5n, Z = S^,. 

From this setup, the following components can be calculated: 

2 

Guu = — {f,u — 2//i), 
r 

= ^ (^'^/.^ + + Afg) , 
2 

Gvv = — {g,v - 'igd), 

r 

Goo = -4^ (du + ^ 



T — —W'^ 
47r 

Tuv = Y^(5), 

rp ^_ y2 

-'-vv — . ) 

2 

Te9 = ^WZ-r^ViS), (3.5) 



where 



v{s) = vm (3.6) 

From the equation of the scalar field, we get the following equation: 

rZ^u + fZ + gW + TTa^rV' {S) = 0. (3.7) 

Note that, V' {S) = dV{S)/dS. 

Finally, we use the Einstein equation, 

Gfj,u = SnTf^,^. (3.8) 

Including this, we can list all equations for our numerical simulations. 
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1. Einstein equations: 



d,u = = ^ + -7^ - WZ, 
g^y = 2dg - rZ'^, 

9,u = f,v = -ll-^ + 27Ta'rV{S), 
r 4r 

t2 



/,„ = 2fh - rW'. (3.9) 



2. Scalar field equations: 



Z,u = W,^ = WV'{S). (3.10) 

r r 

3.2 The initial value problem and integration schemes 

We prepare a false vacuum bubble along the initial ingoing surface, where the outside is 
flat background. One can interpret that a combination of fields is generated in the almost 
flat background via quantum tunneling. We need initial conditions for each function on 
initial n = = and v = v\ = Q surfaces. There are two kinds of information: geometry 
(q, r, f, h, d) and matter (S", VF, Z). 

On the geometry side, we have gauge freedom to choose the initial r function; although 
all constant u and v lines are null, there remains freedom to choose the distances between 
null lines. We choose r{u,v\) = ur^Q + tq and r{u\,v) = vr^Q + tq. Here, we fix tq = 10. 
Then, g{u\,v) = r^Q and f{u,v\) = r^o are naturally obtained. We assume that the 
asymptotic outside is flat: a{u\^v-i) = 1. Since the mass function(M(u, v) = (r/2)(l + 
Ar^uf^y/a^)) should vanish at the initial surface (tti,u), we can choose r„o = ~l/2, 
Vvo = 1/2. 

On the matter side, we fix S{u\^v) = and S{u,Vi) will be defined in the next sub- 
section. Then, one can calculate S{u,Vi), W{u,Vi), S{ui,v), and Z{ui,v) for initial states. 
Then, as one fixes S{u,Vi), from the Einstein equations, one can obtain a{u,Vi) from 
2fh = rW^ ( since r-^^ — at the initial surface). And then, the other functions can be 
evolved using equations on a^uv, f,uu or r^^, and S^uv 

We can choose two integration schemes. First, we can get a from d, r from the equation 

for 

'^,vv-) S'lid. S from Z. SGCond, wg Ccin get oi from t from tliG GQUBtion for Ty^y^^ and 
S from W. We call the former u-scheme, whereas the latter u-scheme. We mainly used 
the t^-scheme. However, these results should be same. We compared them to check the 
consistency of simulations in Appendix Here, we used the 2nd order Runge-Kutta 



method |25| 



3.3 Free parameters 

Now we specify initial parameters. 

First, we specify the potential function. We want a potential that has two stable 
minima, where the true vacuum of a field value $ = has vacuum energy and the false 
vacuum of a field value $ = $0 has vacuum energy A. 
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Second, we specify the initial field configuration Uj). For simplicity, we assume 
that the inner part has a field value <I>o and the outer part has a field value 0, i.e. the true 
vacuum. However, to go beyond the thin shell approximation, we need a transition region 
from the true vacuum to the false vacuum. So, we define the function ^{u,Vi) as follows: 



U < -UshcU, 

$oG'(n) Ushell < -U < Ushell + At(, (3.11) 
^0 ^ishell + Au<U, 



where G{u) is a pasting function which goes from to 1 by a smooth way. We choose G{u) 

by 



G{u) = sin 



7r{u - Ushell) 

2Au 



(3.12) 



and we choose Usheii = 5 for all simulations in this paper. ^ 

Now, we can specify all free parameters of our simulation. 

1. A potential V{^) 

2. Thickness of the transition region Au 

3. A field value of the false vacuum <I>o 



^Here, the initial energy of the shell is proportional to (<l>o/Au)^ since the energy of shell depends on the 
gradient of the field, as we can see in the Tuu component. Note that if the radial function of the transition 
region increases, we interpret that the shell expands; if the radial function decreases, we interpret that the 
shell collapses. Since < and r^v > at the outside of the shell, if the transition region approaches the 
out-going null direction, the shell expands; if the transition region approaches the in-going null direction, 
the shell collapses. 
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( ^11 = 0.1, So = 0.1, A = 0.001, Vpd>) 
r 




10 20 30 40 50 




Figure 6: Simulations of r and S for (Au = 0.1, S'o = 0.1, A = 0.001) with potential Vpoiy(5'). 
Black curves of upper diagram are contours of r. The lower diagram is for S. Here, the yellow 
region is in the false vacuum, and the blue region is in the true vacuum. This shows the collapsing 
shell solution. 



4. Vacuum energy of inside A 

From now, for convenience, we use S = -v/47r<I> rather than ^ for numerical simulations; 
but, of course, we can change both conventions easily. Therefore, if one fixes a potential 
V{S), then three parameters {Au,So = \/A^^o,A) fully defines a simulation. 

4. Causal structures and physical issues 

In this section, we will classify 4 types of solutions. First, as we discussed in the thin shell 
approximation in Section ^, we observe a collapsing shell solution. We call this Type 1. 
Second, by giving sufficient energy to a shell, we will get an expanding shell solution. We 
call this Type 2. However, our simulations do not contradict with the expectations of the 
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Figure 7: Variation of the thickness of the shell Au. As Au decreases, the shell expands, and field 
values of the inside false vacuum become unstable. 



thin shell approximation, since field values of the inside false vacuum become unstable. 
To induce inflation, we will prepare N matter shells and exotic matter shells. In this 
setup, an inflating shell and creation of a bubble universe are possible. We call this Type 3. 
Finally, by tuning the initial parameters, we observe a transition from Type 2 to Type 3; 
between these two types, we could flnd another solution which we call Type 4. 

4.1 Type 1: a collapsing shell 
4.1.1 Collapsing shell solutions 

We begin with the following potential Vpoiy(^) of polynomial form: 



(4.1) 



This potential has a stable minimum around $o with vacuum energy A. The three param- 
eters A, ^>o, and A define one specific potential. For example, if one chooses A = 10^, Sq = 
0.1, and A = 0.001, then the following potential V5poiy(5') is obtained (Figure ^). As we 
change So and A, by tuning A, we could tune the ratio between the height of the unstable 
equilibrium and the height of the false vacuum to be ~ 4.5. 

First, we calculate (An = 0.1, Sq = 0.1, A = 0.001). Figure ^ shows the results. 

The upper diagram of Figure |6| is for the function r. Gradients of each contour lines 
are changed around u = 5, since there is the transition region from the true vacuum to the 
false vacuum. Though the gradients are changed from lower to upper region, each contour 
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Figure 8: Variation of the field value of the inside false vacuum 5*0. As Sq increases, the shell 
expands, and field values of the inside false vacuum become unstable. 



line in the upper region is almost straight and parallel. This implies that the upper region 
of u = 5 does not inflate.^ 

The lower diagram of Figure |6| is for the function S. The yellow region is the false 
vacuum, and the black region is the true vacuum. One may notice that there are small 
fluctuations in the yellow region and the black region. However, even though there are 
fluctuations, the field value of the inside false vacuum region is almost constant. Thus, the 
false vacuum is quite stable. 

One can notice that the false vacuum bubble collapses along a time-like direction. 
Then, eventually it will form a black hole, though we cannot see the black hole in our 
simulation. This result is demonstrated in Figure We call this solution Type 1. This 
result is consistent with dSA — Schn solution of the thin shell approximation. 

4.1.2 Transition from stable to unstable field values 

As we tune the three parameters (An, 5*0 = VAtt^q^ ■^): we can observe the change of the 
collapsing shell solution. Note that all r diagrams or causal structures are similar, and 
hence we omit r functions and present only diagrams of the field S. 

^Here, we could not simulate beyond r — 0, since it gives a mathematical singularity of our equations; 
and hence, from r = point, there is a cutoff line where its beyond is impossible to calculate. This problem 
comes from the 4-dimensionaI spherical symmetry. However, this should be resolved as we restore full 
4-dimensionaI gravity and hence it may not be a fundamental limitation. 
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Figure 9: Variation of tiie vacuum energy of tlie inside false vacuum A. As A decreases, the shell 
expands, and field values of the inside false vacuum become unstable. 



First, we have varied the thickness of the sheh Au (results shown in Figure |^ and 



Figure ^): 








• 


(An = 


0.5, So = 


0.1, A = 


0.001), 


• 


(An = 


0.1, 5o = 


0.1, A = 


0.001), 


• 


(An = 


0.05, So -- 


= 0.1, A = 


= 0.001), 


• 


(An = 


0.025, 5*0 


= 0.1, A 


= 0.001) 


• 


(An = 


0.01, 5o = 


= 0.1, A = 


= 0.001). 



Note that the energy of a shell is approximately proportional to Sq/Au^ since the energy- 
momentum tensor contributes order W^. Therefore, small An and large 5*0 implies more 
energetic shells. As the thickness of the shell becomes smaller and smaller, the shell tends to 
collapse more slowly. In (An = 0.05, So = 0.1, A = 0.001), fluctuations outside of the shell 
(i.e. near the true vacuum) becomes larger than fluctuations of (An = 0.1,5*0 = 0.1, A = 
0.001). In (An = 0.025, 5o = 0.1, A = 0.001) and (An = 0.01,^0 = 0.1, A = 0.001), the 
shell seems to be expanded and slowly rolls down to the true vacuum. Note that the inside 
region (high u) rolls down faster than the region near the shell (especially visible for the 
case An = 0.025). 

Second, we have varied the field value So of the inside false vacuum region (Figure |8|): 
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Figure 10: Transition from Type 1 to Type 2, as the energy of the shell decreases. 



• 


(An = 


0.1, So 


= 0.08, A 


= 0.001) 


• 


(An = 


0.1, So 


= 0.1, A = 


= 0.001), 


• 


(Au = 


0.1, So 


= 0.15, A 


= 0.001) 


• 


(An = 


0.1, So 


= 0.3, A = 


= 0.001). 



As the field value becomes larger and larger, the shell slowly collapses. In (An = 0.1, 5*0 = 
0.15, A = 0.001), the shell expands, but since the inside field values are unstable, the 
field quickly rolls down to the true vacuum and oscillate around S = 0. The case {Au = 
0.1,^0 = 0.3, A = 0.001) is again an expanding solution with slow-rolling field values. 

Third, we have varied the vacuum energy A of the inside false vacuum region (Figure ^ 
and Figure ^) : 



• (An = 


0.1, 


So 


= O.L 


A 


= 0.05), 


• (An = 


0.1, 


So 


= O.L 


A 


= 0.01), 


• (An = 


0.1, 


So 


= O.L 


A 


= 0.001) 
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• (An = 0.1, = 0.1, A = 0.0001), 

• (An = 0.1, = 0.1, A = 0.00001). 

When there is large vacuum energy, the shell tends to collapse quickly. As one has small 
A, the shell tends to expand (An = 0.1, 5o = 0.1, A = 0.0001) and the inside field value 
eventually becomes unstable (An = 0.1,5*0 = 0.1, A = 0.00001). According to the thin 
shell approximation, as the vacuum energy decreases, the critical energy which determines 
collapse or expansion becomes decreases Therefore, as A decreases, the expanding 
solution will be easily obtained and hence our observation is consistent in terms of the thin 
shell approximation. 

Therefore, one can conclude that as the energy of a shell increases, the shell tends to 
expand. This result is consistent with the arguments of the thin shell approximation. We 
observed basic tendencies as follows. First, if the shell has sufficiently low energy than a 
critical value, i.e. large An, small Sq, and large A in a certain limit, then the shell tends 
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0.15 0.2 



Figure 12: An expanding shell with unsta- Figure 13: Kraii(S') with {wi — 0.05, W2 = 
ble field values (Type 2). 0.025, Abump = 0, 0.0035, 0.01). 

to collapse. In other cases, the shell will tend to expand. Second, if the shell collapses, 
the inside false vacuum region is stable; while if the shell expands, the inside false vacuum 
region is unstable and the field values slowly roll down to the true vacuum. In the latter 
case, the rolling begins from the inside to the outside. We include a transition diagram 
between two extreme cases (Figure p^). 

4.2 Type 2: an expanding shell with unstable field values 

4.2.1 Expanding shells with unstable field values 

From the discussions of the previous section, we observed that as the energy of the shell 
increases, the shell tends to expand. In this section, we observe details of the expanding 
unstable field value solution. We will call this Type 2. 

We simulate the following parameters: (An = 0.1,5*0 = 0.3, A = 0.001) (Figure ^). 

The upper diagram is a contour diagram of the function r. The causal structure is 
similar to Type 1. The lower diagram is of the function S. The field value slowly rolls 
down to the true vacuum. The inside region (high u) rolls down faster than the outside 
region. This result is schematically shown in Figure |l^. 

4.2.2 Dependence on potentials 

One question is whether the unstable behavior of inside field values depends on the potential 
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structure or not. For this purpose, we introduce the fohowing form of potential (Figure 13): 



^wall(S) 



r 

A 
2 



1 



COS 



A + 
A 



SVr 

Wl 

1 — COS 



W2 



s <0, 

0< S <wi, 

Wl < S < Wl + 2W2, 

Wl + 2w2 < S. 



(4.2) 



Here, we have used wi = 0.05 and W2 = 0.025. Therefore, as we increase Abump, the 
barrier of potential becomes higher and higher. If the field dynamics is not changed via 
these different heights of barriers, it will imply that the dynamics of unstable field values 
do not sensitively depend on the specific potential structure. 



We calculated the following cases (Figure 14): 



Abump — 0, 

Abump = 0.0035, 



A 



bump 



0.01, 



where the other conditions are fixed by (An = 0.1, 5o = 0.3, A = 0.001). We observed 
that their field dynamics are almost the same. Therefore, the unstable behavior seems to 
depend on just the properties of the shell and does not sensitively depend on potential 
structures. 
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4.2.3 Stability analysis 

One interesting question for an expanding bubble is whether there is inflation or not. In 
Type 2, there is no inflation. Then, what is the reason of this? 

In this paper, we will focus on the case when the fleld value moves slowly. Note that it 
does not mean the same as conventional slow-roll inflation, it just means that the change 
of the fleld amplitude via dynamics of the field is sufficiently small. ^ To induce infiation, 
we need to maintain the fleld values of the inside false vacuum, and as we observed, it does 
not sensitively depend on the form of the potential. Therefore, for this analysis, we use a 
simplifled potential as shown in Figure 15. Here, $,n|sheii — Sq/Au and V^^^n = A/5o. 

We want to see whether the near of the shell induces 
inflation or not. We will discuss brief outline of this section. 
First, to induce the stable field values, the back-reaction 
of the fleld should be sufficiently small, i.e. 



V(S) 



A 







< 1 



(4.3) 



So 



S(u) 



So 



Au 



Figure 15: A simplified poten- 
tial and initial conditions. 



should hold, where An is the thickness of the shell and Av 
is the required time to observe infiation. Second, to observe 
the effect of the vacuum energy, we require a ~ 1; if it 
does not hold, then all contributions of the vacuum energy 
will be suppressed by small a. These two constraints will 
give sufficient forms of the initial conditions Au,^o, and 
A. Finally, we will see that, even though we require the 
conditions, it is difficult to obtain infiation, since we have 
to require sufficiently long Av. 

Let us observe details. If one integrate the scalar field 
equation along du, it contributes only Au. Then the initial 



velocity of field becomes, approximately 



JQ 

r 



2 A ^ 

vra ttAu 
•^0 



(4.4) 



around the shell. Therefore, if the field value of the false vacuum or the vacuum energy 
of the inside false vacuum is too large, the inside field values are unstable. To make the 
field values stable, we need small 5*0 and small Au limit with A An ^ S'o ^ 1; we need to 
choose a proper relation between 5o, An, and A. 

Since r^^u = and r^^ = — 1/2 along the initial constant v line, 

2 



-(In a) 



holds. Then approximately. 



a 



exp 



An 



q2 

._fo_ 
An 







(4.5) 



(4.6) 



Via the instability of fast-rolling fields (see Appendix A.l), calculations of fast-rolling fields require 
too large computing power, and the authors could not obtain any evidence whether there exist fast-rolling 
inflation in our setup. 
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Therefore, one may need the condition 

^shell|| < 1. (4.7) 

If it does not hold, then all terms which depend on potential V are decoupled from all 
equations. If the field amplitude Sq is sufficiently small, one can assume the condition. 
Now let us assume a ~ 1. Note that, this condition implies that the curvature around 
the shell is sufficiently regular; if a ^ 1, then the curvature i? ~ 1/a^ cannot be regular 
around the shell. 

If one compares two values 



IS'ol kshelll' 

this ratio will be sufficiently small as one assumes A-y/rgheii ^ 1 and r,^ ~ 1. Here, Av is 
the required time for an observation of inflation. Also, we need to compare two values 

|g2| ~ |Ar-sheuAi;|. (4.9) 

Here, A ~ lAgheii ^ natural guess, and Au/rgheii <C 1 is needed. Then, |$^uj;At(Ai;|/|$| 
becomes sufficiently small and we can justify a stable field values. 

Therefore, we can suggest some necessary conditions for a stable vacuum: 

1. a ~ 1, or rsheii>S'o/AM < 1, i.e. the curvature around the shell is sufficiently regular, 

2. r ^ ~ 1, 

3- A ~ l/r,\,„ 

4. Aw/rshcu <. 1. 
Then it is natural to assume the following initial conditions: 

An~e^, So — , A ~ (4.10) 

with large rgheii and small e limit. 

However, this situation is difficult to implement in real situations. If we write the 
equation for r^uv and use the necessary conditions, one can obtain the following equation 
near the inside of the shell: 

r„. = /,. = -^-^ + 27raVy 
f 1 

— + 27rrA. (4.11) 

r 4r 

If there is inflation, then initial / = —0.5 will increase to 0. As one integrate both sides 
along V, the first and the second term are almost same order before the beginning of 
inflation. Therefore, more important contribution of / comes from the third term. And, 
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near the very thin sheh, it is reasonable to choose r ~ Cv + rsheib where C is a constant. 
Then approximately, 



f{v) - f{vO ~ C (-^) ' + (^) . (4.12) 

In this limit, Av for f{Av) = is on the order of Tsheii- Therefore, Au/rsheii may not be 
sufficiently small. This is a basic intuitive reason why in our setup it is difficult to induce 
inflation. 

4.3 Type 3: an inflating shell 
4.3.1 Physics of exotic matters 

Before we discuss a new type of solutions, we comment on the physics of exotic matters. 
If a matter violates the null energy condition, it is called a phantom matter or an exotic 
matter. This kind of matter were discussed in many contexts including general relativity 



| 2q| and cosmology [pTj] . 

We prepare the following Lagrangian for an exotic matter field ^: 

C = +^.a^.bg''' + '2V{^). (4.13) 

From this Lagrangian, we can derive the equation of motion for the scalar field: 

*;a6ff"'-^'(*) =0. (4.14) 

Also, the energy-momentum tensor becomes 

Tab = -*;a^';fe + ^9afe (^;c^';d/'' + 2V {^)) . (4.15) 

Note that the potential V{^) contributes to the negative vacuum energy. Since our aim is 
to induce inflation, for convenience, we choose the potential by V{^) = 0. 

According to previous researches, if there exist exotic matter, a static wormhole, a 
warp drive or a time machine may be possible ||2^. However, these analyses were based on 
a pre-existing metric ansatz which may not be justified from an almost flat background. In 
our setup, we do not assume a strange geometry from the initial condition; we begin with 
an almost flat background and assume a combination of flelds. 

By assuming the violation of the null energy condition, an expanding and inflating 



bubble solution may be justified without Farhi-Guth-Guven tunneling |12|. However, still 
the dynamical causal structure of these inflating solutions are not known. In the following 
sections, we will discuss the correct causal structure of the expanding and inflating bubble 
solution. Here, a wormhole is dynamically generated along the shell. 

4.3.2 The A^-shell bubble 

To induce inflation, one may need to hold the field values for sufficiently long time, i.e. we 
should control the rolling of the inside field values. We need some conditions: (1) the force 
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term V'{S) should be suppressed; (2) the vacuum energy of the inside of the shell should 
be maintained; (3) the curvature around the shell should be sufficiently regular. 

We introduce a brief outline of this section. One observation is that if the amplitude 
of a field is sufficiently small, then we can expect that the contribution on V'{S) can be 
sufficiently small. However, the corresponding vacuum energy also becomes small. To 
solve this problem, we introduce a number N of shells. Then we can dilute the force term 
maintaining the vacuum energy as a constant value. However, in that case, the curvature 
around the shells cannot be small, since the contribution on the energy-momentum tensor 
becomes order s/N. To regularize the curvature around the shells, we will introduce a 
number of exotic matter shells. Then we can induce sufficient setup to induce inflation. 

We discuss the details as the following. Let us assume N scalar fields 0^ with potential 



N 



VpM4{(t>i) = ^\J2<t>fj, (4.16) 

where (f>i arc scalar fields and A is a constant. If all fields are coherent and have the same 
amplitude, then V^hi4(0i) = NX(f)j holds. Here, the field equations are 

('/'0;afeff"'' - 4A0f = 0, (4.17) 

where <pi ~ with N fields. Also, assume N exotic matter fields Q. Then, the 

equations of motion are 

(0);a65"' = 0. (4.18) 

Again, we assume Cj - l/iVV4, and assume that all (j are coherent with the same ampli- 
tude. 

Let us assume initial conditions for (^j and by 



(pi{u, Vi) 



U < 7ishcll, 

(</.o/iVi/4)G'(n) Usheii < < «sheii + Au, (4.19) 



U < ■Ushell, 

Cjiu, Vi) = { if3(t>o/Ny^)G{u) Wshcu <u< Ushell + Au, (4.20) 
((/)o/iVV4 ^^^^^^ + Au<u. 



Here, we define 1 - /J^ = l/^/N. 
Then approximately. 



' (4.21) 



l(A);.i9°''l AT'/" 

with large A^. Thus, the force term can be sufficiently suppressed in this setup. 
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Contributions to Einstein equations should be carefully checked. 

2 

V T OL 

(lna),„^ = + 4^ - ^'^^{<^i)Mi).v + 47riV(C^.)^„((^.)_„ 



r = 2r - 47rr(7V(</,i)2^ - 7V(0)'J = 2r,„^ - 47rrViV$2^ + 47rrViV$'5„ 

n« = 2r,„^ - 47rr(iV(</.i)2^ - iV(C,)fj = 2r,„^ - 47rrViV$f„ + 47rrViV$'^„, 



2 2 

r uv = — ^ — — I- 27rQ; rNXcp, = — ■ — ■ — h zira r A$ , 



r r JN 2 



^, _^_r^_ (4.22) 



where we define effective fields <^{u,v) = N^/^(j)i{u,v) and (^'{u,v) = N^f^C,j{u,v). We 
assume the following initial conditions: 

U< Ushell, 

^{U, Vi) = { 4>oG{u) lishcU <U< «shell + Au, (4.23) 
UsheW + Au<U, 



Now they are equivalent to the following scheme: 



U < ■UshcU, 

P(poG{u) Ushell <U< Ushell + An, 
P(pO Mshell + An < U. 



(4.24) 





= hv 




9,v 




~ f,v 






- 






= K 



2dg - T\fNZ^ + r f 1 \ 

— + 27ra r-^S , 

r 4r Sq 



r r 
/Z' gW 



1 A 



(4.25) 



u< nshcu, 

5, S'{u, Vi) = { SoG{u) UshcU < U < Ushell + An, 

So nsheii + An < u, 



(4.26) 
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{/Ill = 0.3, So = 0.1, A = 0.01, N=10«, VpM) 




Ig 20 30 40. ,5ft 



Figure 16: r diagrams for {Au = 0.3, So = 0.1, A = 0.01, iV = 10^) and (Au = 0.3, So = 0.1, A = 
0.05, N = 10^) with Vp\ii4{S) potential. White curves are r^u — horizons. Blue region in the lower 
diagram is space-like future infinity. This shows an inflating shell solution. 



since the equation for is linear and hence one can re-scale about factor p. Here, So = 
\/47r$o is the field amplitude and A = (47r)^A/S'Q is the vacuum energy of the inside of the 
shell. 

Let us observe a naive expectation of this formulation. The field equation for $ is 
effectively free via 1/y/N term. Therefore, effectively, $ ~ ^frec + 0{1/^/N) is obtained. 
Contributions to the energy-momentum tensor is on the order of 

i.e. the leading term contributes to the energy-momentum tensor by order ^/N, where 
the second term contributes order 1 which comes from the V'{S) effect. However, the 
first contribution will be canceled by A'^ exotic matter shells, and hence their contributions 
become order 1. Therefore, even if we assume N ^ oo limit, we cannot naively assume $ 
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Figure 17: An inflating sliell solution (Type 3). After push a terminal matter, the evolution of the 
shell is ended via a black hole, (a), (b), and (c) are schematic diagrams of area for each space-like 
section. 

as a free field, since the second correction term is comparable to the free field contributions. 
If $ is a free field, as long as the initial surface maintains the null energy condition, the 
null energy condition will hold in all regions. However, in our cases, even if we prepare 
the initial energy-momentum tensor to hold the null energy condition, the effect of V'{S) 
is not negligible; as time goes on, V'{S) term may affect to roll down the field, and hence 
eventually the null energy condition will be violated during inflation. 

4.3.3 Simulations of A^-shell bubbles 

Here we observe simulations with the following conditions: {Au = 0.3, 5o = 0.1, A = 
0.01, iV = 10^) and (An = 0.3, 5o = 0.1, A = 0.05, = 10*^) with Vpu^iS) potential 
(Figure [l^ ). In the lower diagram, the blue region (upper right corner) is beyond the 
calculation ability, since the radial function r becomes exponentially large. We can interpret 
this as a space-like future infinity. Therefore, this is an important evidence that we induce 
inflation. We will call this class of solutions Type 3. 

One can observe the r^^ = horizon: an ingoing observer sees increase of area. The 
inner horizon has two parts: one is almost parallel to the ingoing null direction, and the 
other is almost parallel to the outgoing null direction. The former corresponds to the cos- 
mological horizon of the de Sitter space. The latter corresponds a region where area begins 
to increase for an ingoing observer: it is similar with the throat of a wormhole. There- 
fore, we can say that, during inflation, a wormhole is dynamically generated (schematically 
shown in Figure [l^). One can observe that the former part is always time-like, while the 
latter part begins a space-like direction and finally becomes a time-like direction (in fact, 
the latter is almost null). 
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{Aii = 0.3, Su = 0.1, A = 0.01, N=IO*, Vphi4) 



1 




i 







Figure 18: Function 5*, S' , T^u, and T„„ of {Au = 0.3, 5*0 = 0.1, A = 0.01, N = 10^) with VpY^aiS) 
potential. Yellow regions of Tuu and Ty^ have field values more than 10"'* and 10"^; skyblue regions 
of Tuu and Ty^ are less than —10"^ and —10"^; the sign of the energy- momentum tensor is changed 
around black region. 



(^11 = 0.3, Sii = 0.1, A = O.Oi, .N=IO^ V|iW4) 
5.4 r 




Figure 19: r„„ and T^v around the shell. 



Note that, in Figure 17, we pushed a terminal matter to induce a black hole. We do 
not need that stage, but to make an end of the structure and discuss the information loss 
problem in the final section, we artificially inserted the process. 

We observed that the wormhole throat can be generated around the mass shell. Hence, 
we do not need to assume tunneling from the outside to inside of a Schwarzschild wormhole. 
In fact, if the null energy condition is violated, this behavior may happen. However, 
its causal structure was not known, since the structure is related to the thick transition 
layer. We observed a dynamical generation of a wormhole around the transition layer, and 
observed the causal structure. 

Also, we observed the field S and S", as well as the energy- momentum tensor T^u and 
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Tyy (Figure 18). Diagrams of S and S' are almost similar, but S rolls down more quickly 
than S'. One can see black bands from T^u and T„^. Note that, if T^u or T^^ is less than 0, 
it implies the violation of the null energy condition. The violation of null energy condition 
seems to begin around the inside of the shell (Figure ^9|). And the exotic matter becomes 
dominant as time goes on. Therefore, we highly suspect that our inflating shell solution 
requires a violation of the null energy condition. 

4.4 Type 4: a bursting shell 

4.4.1 Transition from Type 3 to Type 2: a bursting shell solution 

As we vary the simulation parameters, we could find another kind of solutions between 
Type 3 and Type 2. 

First, we have changed vacuum energy of the inside A: 



• (An = 


0.3, ^0 


= 0.1 


A 


= 0.05, iV = 


io«), 


• (An = 


0.3, ^0 


= 0.1 


A 


= 0.01, iV = 


io«), 


• (An = 


0.3, 


= 0.1 


A 


= 0.005, = 


= io«), 


• (An = 


0.3, So 


= 0.1 


A 


= 0.001, = 


= 10^). 



Second, we have changed the number of shells A^: 
• (An = 0.3,50 = 0.1, A = 0.01, A^ = 10^), 
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Figure 21: A bursting sliell solution (Type 4). After push a terminal matter, the evolution of the 
shell is ended via a black hole, and a Cauchy horizon is generated, (a), (b), and (c) are schematic 
diagrams of area for each section. 

• (An = 0.3,5'o = 0.1,A = 0.01,iV = 10*^), 

• (An = 0.3,5o = 0.1,A = 0.01,iV = 10^). 

Figure m shows the transition from Type 3 (An = 0.3, S'o = 0.1, A = 0.01, iV = 10*) 
to Type 2 (An = 0.3,50 = 0.1, A = 0.001, iV = 10*). Between these two limits, we could 
find interesting structures: we will call these Type 4 (Figure ^l|) . 

The cases (An = 0.3,50 = 0.1, A = 0.01, iV = 10*^) and (An = 0.3, 5o = 0.1, A = 
0.01, N = 10^) show the behavior when N decreases. Inflation is suppressed, but the loca- 
tion where inflation begins is invariant. However, in (An = 0.3, S'o = 0.1, A = 0.005, A/^ = 
10*), we observe small A limit, and the beginning of inflation is shifted. 

If one compares stability of field S (Figure ^2|), as N decreases, inside field values 
become unstable and roll down to the true vacuum and roll around 5 = 0: (An = 0.3, 5o = 
0.1, A = 0.01, = lO''). 

Figure ^ shows Tuu components. Yellow regions are greater than 10~^ and skyblue 
regions are less than — 10~^. In the case of Type 3 solutions (e.g., (An = 0.3, 5o = 0.1, A = 
0.05, = 10*) or (An = 0.3, 5o = 0.1, A = 0.01, A^ = 10*)), violating region of the 
null energy condition is quite wide. However, in Type 4 solutions (e.g., (An = 0.3, 5o = 
0.1, A = 0.01, A^ = 10^), (An = 0.3, 5o = 0.1, A = 0.01, A^ = 10^), or (An = 0.3, 5o = 
0.1, A = 0.005, A^ = 10*)), Tuu is almost globally positive except near the shell. Therefore, 
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Figure 22: Field S for some conditions. As N decreases, the inside field more quickly rolls down. 



the violation of the null energy condition seems to be essential, but to see r „ = horizons. 



the violation is needed just a narrow region around the shell (Figure 24). 



(^u = 0.3, So = 0.1, A = 0,005. N = IO», Vphii) 




Figure shows T^y components. The 
beginning of inflation seems not to be related 

t""™ ^^^^^^^^^^^^^^^^HBHI^H to the sign of T^^ and may not be important 
OIWO ^^^^ 
-a.oiso for inflation. 

Finally, we remark on the behavior of 
r^u = horizons. For the Type 3 case, the 
inner r^u = horizon (the analogous horizon 
of the de Sitter space) is always time-like but 
for the Type 4 case, the inner r ^ — horizon 
bends from a time-like direction to a space- 
like direction. Also, the outer r^^ = horizon (the throat of a wormhole) bends from 
a space-like direction to a time- like direction (Figure ^). These changes come from the 
properties of r^uv = f,v around the horizon. Since r^u = 0, it becomes 

2 

r „,J. „=o = + 2WrV{S). (4.28) 
4r 



Figure 24: T^u around the shell for (Aw 
0.3, So = 0.1, A = 0.005, N = IQ^). 



If r^uv > 0, then the outer r^u = horizon is space-like and the inner r ^ — horizon is 
time-like; if r^„^ < 0, then the outer r^^ = horizon is time-like and the inner r^u = 
horizon is space-like. The sign of r^^v is positive or negative if and only if 

87rrV(5) - 1 (4.29) 
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is positive or negative. Therefore, if field values are sufficiently large and the vacuum energy 
is sufficiently large, the outer = horizon is space-like and the inner r — horizon 
is time-like. However, as time goes on, the field values slowly roll down and then if the 
vacuum energy becomes sufficiently smaller than a critical value, the outer r^u = horizon 
becomes time-like and the inner — horizon becomes space-like. In other words, if 
the field values quickly roll down to the true vacuum, the horizons tend to disappear and 
inflation ends. This is the basic physical difference between Type 3 and Type 4. 



These changes are schematically shown in Figure 26. 



5. Discussion 



We considered a false vacuum bubble inside of an almost fiat background. Here, we as- 
sumed that the size of the bubble was smaller than the background horizon size, and of 
similar order as its own horizon size. A traditional approach on this problem is the thin 
shell approximation. In this paper, we extended our methods to beyond the thin shell 
approximation . 

To summarize our discussion, firstly we will focus on the meaning of beyond the thin 
shell approximation. Secondly we will discuss about some speculations of a bubble universe 
and unitarity issues. 
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( 1= 0.3, So " 0.1, A = 0.05, N^IO", VphM) (Zlu = 0.3. .So = 0.1. A= 0.01. N=10f>. V,,l.ii) 




Figure 25: T^^, plots for some conditions. Yellow regions are greater than the upper bound, and 
sky blue regions are less than the lower bound. 

5.1 Beyond the thin shell approximation 

We discussed previous results of the thin shell approximation. The thin shell approximation 
assumes the inside false vacuum region as the de Sitter space and the outside true vacuum 
region as the Schwarzschild space. According to the thin shell approximation and the null 
energy condition, if a shell is outside of a Schwarzschild black hole it cannot evolve to a 
bubble universe unless one introduces the Farhi-Guth-Guven tunneling; the only allowed 
solution is dSA — SchD- 

We extend the analysis beyond the thin shell approximation by using numerical meth- 
ods. Two essential points of beyond thin shell approach are the added field dynamics and 
the thick transition layer. If a shell has sufficiently low energy, as expected from the thin 
shell approximation, it will collapse (Type 1). However, if the shell has sufficiently large 
energy, it tends to expand. However, via the field dynamics, the inside vacuum slowly rolls 
down to the true vacuum (Type 2). Moreover, if we add sufficient exotic matters to reg- 
ularize curvatures around the shell, inflation may be possible without Farhi-Guth-Guven 
tunneling and a wormhole can be dynamically generated (Type 3). By tuning parameters, 
we could find transitions between Type 1 and Type 2, as well as between Type 2 and 
Type 3. Between Type 2 and Type 3, we also find another class of solutions (Type 4). 

We can conclude our new contributions on this issue. First, finding of Type 2 is a 
new result, thanks to going beyond the thin shell approximation. Second, we discussed the 
physical reason why inflation is difficult to induce without the violation of the null energy 
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Figure 26: Transition from Type 3 to Type 2. 



condition. Third, the exact causal structure of Type 3 is also a new feature. Some authors 
know that the violation of the null energy condition may allow an expanding inflating 
bubble, but the causal structure was not known exactly since the structure is related to 
the thick transition layer. We observed the causal structure of an inflating bubble as 
the null energy condition is violated, and it accompanies the dynamical generation of a 
wormhole. Fourth, finding of Type 4 is also a new result. Here, we observed the properties 
of r^u = horizons and the energy-momentum tensor; properties of = horizons are 
related to the potential V, and to induce inflation, T^u < seems to be the sufficient 
condition. Fifth, we could observe a continuous change of each types as one tunes initial 
parameters. In fact, the important parameters are just on shells: An, ^q, A, and effects of 
the other parameters may not be important. 
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5.2 Generation of a bubble universe: discussion on the information loss problem 

If the false vacuum bubble is greater than the background horizon size, even though the 
bubble is separated from the background, it is not so meaningful since an observer of a 
scattering experiment is always inside of the bubble. However, in this paper, we discussed 
a bubble which is sufficiently smaller than the background horizon size. Therefore, in prin- 
ciple, it can be discussed in the context of scattering experiments, and thus the generation 
of a bubble universe should be discussed in the context of the information loss problem 

The Type 3 and Type 4 solutions induce a separation between the outside true vacuum 
region and the inside region. The inside will have a second asymptotic region; Type 3 case 
is quite clear in this point. Then, information loss is inevitable unless duplication of 
information happens.^ If the background is an anti de Sitter space, one may notice the 
AdS/CFT correspondence that implies unitarity In this context, some authors said 
that Farhi-Guth-Guven tunneling should be excluded In the same sense, we can say 
that some assumptions of our setup is inconsistent with unitarity and AdS/CFT. What are 
the exotic assumptions of our setup? Firstly, we assumed the existence of exotic matter 
fields. Secondly, we assumed some special initial conditions of scalar fields and exotic fields. 
Thirdly, a large number of N shells was used to maintain field values. 

We can comment on the second assumption. If a bubble is generated in an almost 
flat background, it can be described by a field combination for the inside false vacuum 
and the outside true vacuum. Here, if the background is a de Sitter space or an anti de 
Sitter space, the tunneling is not excluded via the violation of the energy conservation |5|. 
We found that initial conditions for each types are not quite different; each types transits 
continuously (Figure ^ and ^). Therefore, if certain initial conditions are excluded by an 
unknown reason, the other conditions should be excluded, too; if a combination for Type 2 
is possible in principle, then it is difficult to find a reason why a similar initial condition 
that generates Type 3 should be excluded. 

We can comment on the large number N. To see the r^u = horizon, we do not need 
such large A^. Also, there are some discussions that string theory seems not exclude a large 
number of fields ||3^ . 

The most strange assumption is the existence of exotic matter fields. Now, we can say 
that if one assumes the existence of the exotic matter fields, the creation of a bubble uni- 
verse, dynamical generation of a wormhole and the violation of unitarity become possible. 
We show a clear contradiction between the existence of a certain combination of exotic 
matter fields and unitarity or AdS/CFT. 

Then, the next question is that, does the nature exclude such combination of exotic 
fields? One of physically reasonable processes of an exotic matter or the violation of energy 
conditions come from Hawking radiation For a black hole case. Hawking radiation 

of negative energy tends to inside of the black hole; however, here we need out-going 
negative energy fiux. Therefore, it is unclear whether Hawking radiation may be helpful 
on this issue. However, if we can control the negative energy flux of Hawking radiation 

*Black hole complementarity may be related to this issue. For further discussions, see [H]|28|. 
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Figure 27: Convergence tests of r for the 
condition (Au 0.1, 5o = 0.1, A = 0.001) 
with potential lpoiy('5'). This shows the sec- 
ond order convergence. 



Figure 28: Consistency tests of r by com- 
paring scheme and ^tt — scheme- 




Figure 29: Convergence tests of S for the Figure 30: Consistency tests of 5 by corn- 
condition (Au = 0.1, S'o = 0.1, A = 0.001) paring S't,-schemc and 5„_scheme- 
with potential V^oiy(5')- This shows the sec- 
ond order convergence. 



to the out-going direction, it will definitely make our first assumption reasonable. These 
problems should be discussed later. 
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Figure 31: Convergence tests of r for the 
condition (Aw = 0.1, 5o = 0.3, A - 0.001) 
with potential Vpoiy{S). This shows the sec- 
ond order convergence. 



Figure 32: Consistency tests of r by com- 
piling 7'j;_schcmc and — scheme- 
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Figure 33: Convergence tests of S for the 
condition (Am = 0.1, S'o = 0.3, A = 0.001) 
with potential Vpo\y{S). This shows the sec- 
ond order convergence. 



Figure 34: Consistency tests of 5 by com- 
paring 5'ij_schcinc and *S'ti — schcmG. 



A. Convergence and consistency tests 

In this appendix, we discuss convergence and consistency tests for our simulations of each 
type. Also, we briefly comment on the instability for fast-rolling fields. Here, we check the 
consistency by comparing the u-scheme and the ti-scheme for certain u = constant surfaces. 
We check the convergence by comparing 1x1,2x2 finer, and 4x4 finer simulations for 
certain u = constant surfaces. We observed u = 7.5, 12.5, 17.5 slices. 

We compared two independent evolutions by using equations for r^uu and r^^- Then 
the equation for 

"^^uv remains a constraint equation. Therefore, we checked the following 
quantity to check the constraint equation: 

I - r^ur,v/r - a'^/Ar + 2TTa'^rV{S) - r^uv\ ^-^ 
I - r^ur,v/r\ + \ - a'^/4r\ + \2TTa'^rV{S)\ + \ - r^uv\ 

We checked the constraint equation for v = 10, 20, 40 slices. 
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Figure 35: Convergence tests of r for 
the condition (Am = 0.3, S'o = 0.1, A = 
0.01, iV = 10^) with potential Vphi4(5'). This 
shows the second order convergence. 



Figure 36: Consistency tests of r by com- 
piling ^u— scheme and ^u— scheme • 
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Figure 37: Convergence tests of S for 
the condition {Au = 0.3, S'o = 0.1, A = 
0.01, iV = 10^) with potential Vphi4(5'). This 
shows the second order convergence. 
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Figure 38: Consistency tests of S by com- 
paring 5'ti_schemo and S'ti— scheme- 



Note that simulations in this paper are based on a numerical code of ||l^. And these 
results could be reproduced consistently from an independent code of [^]. 

A.l Type 1 

We checked consistency and convergence for the condition {Au = 0.1, Sq = 0.1, A = 0.001) 



with potential Vpoiy(5'). Figure 27 and Figure 28 show that errors for r is less than 10 



-2c 



Figure and Figure ^ show that errors for S is almost less than 10~^%, but there are 
some peaks which are order of few percents. 

This phenomenon comes from the fast-rolling of the field around the true vacuum. 
Near the true vacuum, as the field rolls around point, a small deviation from may 
look like a big error since we need to divide by 0. This is the instability of fast-rolling 
fields. If its back-reaction to r is sufficiently small, the instability is negligible. Type 1 
holds this case. However, if there exists inflation via fast-rolling fields, one cannot ignore 
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Figure 39: The constraint equation for the 
condition {Au = 0.1, 5o = 0.1, A = 0.001) 
with potential Vpoiy(S'). 



Figure 40: The constraint equation for the 
condition (Au = 0.1,5*0 = 0.3, A = 0.001) 
with potential Vpo\y{S). 



its back-reaction, and to maintain sufficiently small error for S, we need more and more 

finer simulations. This is a practical reason why we consider slow-rolling inflation only. 

Figure 27 and Figure ^ show | (1 x 1) — 

(2x2)1/1(1x1)1 and4|(2x2)-(4x4)|/|(2x 

2)1 for functions r and S, where (n x n) = 

means an n x n finer simulation. In these ' 

figures, for each u slice, two curves are quite | ' 

I 

close and seem to be degenerated. These ' 

show that our simulations converge on the ' 

second order. ' 



Figure 39 shows the constraint equa- 
tion for some slices. The violation of the 
constraint equation is less than 1% for al- 
most all integrated domains (0 < n < 18). 

A. 2 Type 2 




Figure 41: The constraint equation for the con- 
dition (Au = 0.3,50 = 0.1, A = 0.01,7V = 10^) 
with potential VphniS). 



We checked consistency and convergence for 

the condition {Au = 0.1,5*0 = 0.3, A = 0.001) with potential Vpoiy{S). Figure ^ and 
Figure |3^ show that errors for r is less than 10" 



70. Figure |3^ and Figure 34 show that 
errors for S is less than 10~'^%. Also, Type 2 solution converges to second order as well. In 
this case, there was no instability of fields, since the field values of the inside false vacuum 
slowly moves. 



Figure 4C shows the constraint equation for some slices. The violation of the constraint 
equation is less than 1% for almost all integrated domains (0 < u < 18). 

A. 3 Type 3 and Type 4 

We checked consistency and convergence for the condition (An = 0.3,5*0 = 0.1, A = 
0.01, A'^ = 10^) with potential l^hi4('S')- Figure |3^ and Figure ^ show that errors for 
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r is less than 10 %. Figure 37 and Figure 38 show that errors for S is less than 10 



70. 



Also, Type 3 solution converges to second order as well. 



Finally, Figure ^ shows the constraint equation for some slices. The violation of the 
constraint equation is less than 0.1% for almost all integrated domains (0 < n < 18). 
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